Phase diffusion of Bose-Einstein condensates in a one-dimensional optical lattice 
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We investigate the phase diffusion of a Bose-Einstein con- 
densate (BEC) confined in the combined potential of a mag- 
netic trap and a one-dimensional optical lattice. We show 
that the phase diffusion of the condensate in the whole op- 
tical lattice is evident and can be measured from the inter- 
ference pattern of expanding subcondensates after the optical 
lattice is switched off. The maximum density of the interfer- 
ence pattern decreases significantly due to the mixing of the 
phase diffusion appearing in different subcondensates. This 
suggests a way to detect experimentally the notable phase 
diffusion of BECs. 



I. INTRODUCTION 

The remarkable experimental realization of Bose- 
Einstein condensation has generated intensive theoretical 
and experimental investigations for weakly interacting 
Bose gases [1,2] . Although much progress has been made, 
one of the challenging problems to be solved in both 
theory and experiment is how to detect the phase dif- 
fusion of a Bose-Einstein condensate (BEC). The study 
of the phase diffusion is of crucial importance because 
the number of particles realized in all BEC experiments 
is finite. Within the formalism of the mean-field theory, 
it is known that the condensate can be well described by 
a macroscopic wave function (or called an order parame- 
ter) with a single phase [1]. The development of the phase 
at low temperature is governed by the Gross-Pitaevskii 
(GP) equation. However, for a finite system and hence 
for trapped Bose gases, although a macroscopic wave 
function describing a BEC still exists, as is now firmly 
established experimentally, the phase of the macroscopic 
wave function will undergo a quantum diffusion process 
as time develops [3] . Because the phase diffusion reflects 
directly the coherent nature of a condensate, the study of 
the phase diffusion of BECs is not only fundamental but 
also important for applications such as the realization of 
an atom laser [4]. 

For a condensate confined in a trap, the development 
of its phase is determined by the chemical potential of 
the system. Due to quantum and thermal fluctuations, 
the phase diffusion is described by the fluctuations of the 
chemical potential. Even at initial state the condensate 
has an exact single phase, the phase will become unpre- 



dictable due to the fluctuations of the condensate. After 
the realization of Bose-Einstein condensation, much the- 
oretical attention has been paid to the phase diffusion of 
BECs [3-13]. Because for a single condensate it is hard to 
observe the effect of phase diffusion based on the density 
distribution of a Bose-condensed gas, in the experiment 
conducted by the JILA group [14] two expanding con- 
densates are used to produce an interference pattern and 
hence the phase diffusion of the condensates can be de- 
tected. It was found that in this experiment there is no 
detectable phase diffusion on the time scale up to 100 ms. 

Recently, the BECs in optical lattices have attracted 
increasing attention (see Refs [15-19] and references 
therein). Trapping atoms in optical lattices can pro- 
vide very effective and powerful means for controlling 
the properties of a BEC, especially the quantum phase 
effects on a macroscopic scale such as the coherence of 
matter waves [17,19]. In the present work, we investigate 
the effect of the phase diffusion of a BEC confined in the 
combined potential of a harmonically magnetic trap and 
a one-dimensional (ID) optical lattice. In the presence of 
the ID optical lattice, there are an array of disk-shaped 
subcondensates. For each subcondensate, i. e. the BEC 
formed in the corresponding potential well induced by 
the optical lattice, there exist phase fluctuations which 
increase very slowly with the development of time. For 
high enough potential-well intensity, the phase fluctua- 
tions for different subcondensates are independent from 
each other. But when the potential-well intensity is not 
high, there will be a correlation for the subcondensates in 
different wells and hence the situation may be quite dif- 
ferent. One expects that the phase diffusion of the BEC 
in the whole optical lattice can display notable new char- 
acters because of the mixing effect of the phase diffusion 
appearing in different subcondensates. 

The paper is organized as follows. Section II is de- 
voted to the phase diffusion of a single 2D (as well as 
quasi-2D) condensate for the temperature far below the 
critical temperature. In Sec. Ill, the phase diffusion 
of an array of disk-shaped subcondensates induced by a 
ID optical lattice is investigated based on the interfer- 
ence pattern resulting from a superposition of expanding 
subcondensates in different potential wells for different 
holding time of the optical lattice. Sec. IV contains a 
discussion and summary of our results. 
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II. PHASE DIFFUSION OF A QUASI-2D 
CONDENSATE 



8 No « N , the average phase of the condensate is then 
given by 



We consider an array of subcondensates formed in the 
combined potential consisting of a harmonically magnetic 
trap and a periodic optical lattice. The combined poten- 
tial is described by [20-22]: 



4>({No),t) = n({N ),T)t/h. 



(4) 



V 



1 22,1 2 / 2 , 2\ 

-mu x x + -mu ± (y + z ) 



sE R sin 2 (^Y^j 



(1) 



where uj x and wj_ are respectively the axial and radial 
frequencies of the magnetic trap. The last term repre- 
sents the potential due to the presence of the optical lat- 
tice, with A denoting the wavelength of the retroreflected 
laser beam, and s being a factor characterizing the in- 
tensity of the optical potential which increases when in- 
creasing the intensity of the laser beam. In addition, 
E R = 2TT 2 h 2 1 (mA 2 ) is the recoil energy of an atom ab- 
sorbing one photon. For the optical lattice created by the 
retroreflected laser beam, d = A/2 is the distance between 
two neighboring wells induced by the optical lattice. 

The optical potential can be approximated as V e /f — 
J2 k \mZ) 2 x {x - kd) 2 with u x = 2^E R /h [20]. As in Ref. 
[20] we assume that for the parameters appearing in the 
potential (1) the condition u> x « uj± « k R T < ui x is 
satisfied. Thus the magnetic trap is a slowly-varying 3D 
potential which confines the atoms harmonically to over- 
all cigar-shaped distribution. In the ID optical lattice 
the atoms are confined to an array of 2D planes. As a re- 
sult the whole condensate consists of an array of quasi-2D 
subcondensates like disks in shelf. Such array of quasi- 
2D subcondensates has been realized in the experiment 
[20]. 

The investigation of the phase diffusion of the above- 
mentioned BEC can be divided into two steps. The first 
step is to consider the phase diffusion of a single disk- 
shaped BEC. For this aim we develop a general approach 
to calculate the probability distribution of the phase for 
a condensate confined in a quasi-2D harmonic trap. For 
the temperature far below the critical temperature T c , 
the condensate can be described by the order parameter 
with a phase factor <fi (t): 



$ (r, t) = $ (r) e 



i<j>(t) 



The differential of the phase takes the form 
d<p = fi(N ,T)dt/h, 



(2) 



(3) 



where [i(N ,T) is the chemical potential of the system. 
We see that the particle-number fluctuations of the con- 
densate yield fluctuations in the chemical potential, and 
hence lead to the phase diffusion of the condensate. As- 
suming that the mean ground state occupation number 
is (No), in the case of the particle- number fluctuations 



The phase diffusion of the condensate can be described 
by considering the phase difference A(f> (t) = <p(t) — 
4>((No) ,t). From Eqs. (3) and (4), it is straightforward 
to obtain a differential equation on A<p (i) : 



dA<j>(t) _ dti{(N ),T) 



dt 



d(N ) 



AN (t) j%. 



(5) 



In the above expression, AN (t) represents the fluctua- 
tions of the ground state occupation number around (No) 
at time t. Different from the situation considered by Jin 
et al [23] where a small time-dependent perturbation is 
used artificially to create selected collective excitations in 
a BEC, the collective excitations considered here are cre- 
ated spontaneously from the condensate due to quantum 
fluctuations. Thus AiVo (t) changes with the develop- 
ment of time due to the creations and annihilations of 
various collective modes. 

For the temperature far below the critical temperature, 
the collective excitations spontaneously created from the 
condensate play a dominant role in the fluctuations of the 
condensate [13]. We therefore pay attention here to the 
phase diffusion contributed by the collective excitations 
of the system. For a pure 2D Bose-condensed gas, the col- 
lective mode is labelled by two parameters {n, I}. Here, 
the parameter I labels the angular momentum of the exci- 
tation. In this case, one has AN (t) = £„;A7V„; (t). As- 
suming that A(f) n i (t) represents the phase diffusion due 
to the collective mode nl, we have 



dA(j> n i(t) dfi((N ),T) 



dt 



d(No) 



AN nl (t) /n. 



(6) 



Thus the overall phase fluctuations of the condensate 
read A(f>(t) = S„ z A0„ z (t). 

In the above equation, |AA„;| can be regarded as the 
mean number of atoms corresponding to the collective 
mode nl. Based on the Bogoliubov theory [24], |A7V„;| is 
given by 



\AN, 



: (| u 2 nl (v)dV + J v 2 nl (r)dv} f nl + J v 2 nl (r)dV, 

(7) 

where f n i — 1/ (exp(/3e n ;) — I) is the number of the 
collective excitations presenting at thermal equilibrium, 
P = l/(ksT) (T is temperature), and e n i = fru> n i is the 
energy spectrum of the collective mode nl. In the above 
equation, the quantity J v 2 t (r)dV describes the effect of 
quantum depletion, which does not vanish even at T = 
. In Eq. (7), u n i(r) and v n i(r) are determined by the 
following coupled equations: 
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V 2 +V ± (r) - /i + 2gn (r) u„i(r) 



+gn (r) v n i(r) = e nl u n i{r), 



■— V 2 +^ (r)-/i + 2 5 n(r) ) r„,(r) 
+5n (r) u ni (r) = -e n[ v n i(r), 



where n (r) and n (r) are the density distributions of 
the Bose gas and condensate, respectively. V± (r) = 
\mLo\(y 2 + z 2 ) is the 2D harmonic potential confining 
the Bose gas in y and z directions (the second term on the 
right hand side of (1) ), while g is the coupling constant 
describing the interatomic interaction in the condensate. 

By a straightforward calculation the excitation fre- 
quency of the collective mode nl is given by 



uj nl = lu ± (2n 2 + 2n \l\ + 2n + \l\) 



1/2 



while |AiV„;| reads 



|AiV ni | = 



(N mg) 1/2 P nl 



2 v ^h v / 2n 2 + 2n\l\ + 2n + |Z| ' 



(9) 



(10) 



where 



Pnl= [ (l-x 2 )x(H nl fdx/ [ x(H nl fdx, (11) 
Jo Jo 

(12) 



H nl (x) = x\ l \ ]T bjx 2j 

3=0 



with 



b 



and 



the 



efficients bj satisfying the recurrence relation bj +1 
(Aj 2 + Aj + Aj\l\ - An 2 -An- An\l\) / (Aj 2 



CO- 

lh = 

4j + Aj\l\ - An 2 - An- An\l\) / (Aj 2 + Aj\l\ + 8j + 
As mentioned above, the collective excitations consid- 
ered here are created spontaneously due to the quantum 
fluctuations of the condensate. They can disappear as 
time develops and hence have a finite longevity. The 
longevity can be estimated based on the time-energy 
uncertainty relation. For the collective mode nl, its 
longevity is r nl = l/u± (2n 2 + 2n \l\ +2n + \l\) 1/2 . Thus 
at time t » r n i, there is a series of i = t/r n i succes- 
sive creations and annihilations of the collective mode 
nl. Assuming that at the time t the number of the col- 
lective modes created (annihilated) in the condensate is 



{iann)- Thus one has 



+ i 



In this sit- 



uation, the desired probability is given by the binomial 
expression [25] 



PiU) 



(13) 



where j — i cre — i an n- For t » T n i, using the Stir- 
ling's formula, the asymptotic form of the probability 
distribution P n i (A<fi n i,t) of the phase diffusion due to 
the collective mode nl takes the form 



P ni (A<t> nh t) 



1 



x nl 



cxp 



(14) 



^ where the dimensionless parameter r„; is given by 



dn((N ),T) \AN nl \r nl 



d(N ) 



(15) 



From the Gaussian distribution given by Eq. (14), the 
phase fluctuations due to the collective mode nl have the 
form: 



V 2 



nl 



t/r n 



(16) 



Assuming that the collective modes with different nl 
arc created and annihilated independently, one obtains 
that the probability distribution of the overall phase dif- 
ference A<j) (t) is still a Gaussian distribution function. 
The probability distribution of A<f> is then given by 



^(A^,t) = yf cx P (-B(A^ 2 ). 



(17) 



Based on the general probability theory, the parameter 
B can be obtained through the relation S 2 4> = 1/2B = 
(5 2 (p) nl . From the formulas (10), (15), (16) and the 

1 /2 

chemical potential u(N ,T) = (N mg/ir) ' uj± for 2D 
harmonic trap, we obtain the overall phase fluctuations 

S 2 , 



as 



5 2 ^ 



m 2 g 2 Lu±t 



T. 



,+4). 47r 2 n 4 
where the coefficient T takes the form 



^ (2n 2 + 2n\l\ + 2n+\l\) 



3/2 



3rd + 2 



(18) 



(19) 



where (3 n i has been given in Eq. (11). From Eq. (18) we 
see that when the longevity of the collective excitations 
is taken into account, the phase fluctuations are propor- 
tional to the time t, rather than t 2 which is obtained when 
the collective excitations are regarded as quite stable. 

As mentioned above, if uo x « iv± « ksT < ui x , the 
condensate confined in each well has a quasi-2D nature. 
Note that for the array of the subcondensates induced 
by the optical lattice with enough depth, the overlap be- 
tween the subcondensates occupying different wells can 
be omitted. This means that the behavior of the con- 
densate in each well has a local property. For each well 
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the coupling constant g is given by g w 2\^2irh 2 a s / (ml x ) 
[26], which is fixed by a s-wave scattering length a s and 
the oscillator length l x — (h/moj x ) ' in the x-direction. 
Thus from Eq. (18) the subcondensate in each well dis- 
plays the phase fluctuations 



2Tuj±t I a. 



(20) 



From the above result, one has 5 2 (j> <~ lo±uj x . We see 
that the confinement induced by the optical lattices has 
the effect of increasing the phase fluctuations. This is a 
natural character by considering the fact that the con- 
finement has the effect of increasing the particle-number 
fluctuations [27]. 



III. EFFECT OF PHASE DIFFUSION ON THE 
INTERFERENCE PATTERN OF EXPANDING 
ARRAY OF SUBCONDENSATES 

We now go to the second step, i. e. to discuss the phase 
diffusion of an array of subcondensates confined in the 
combined potential of the harmonically trapping poten- 
tial and ID optical lattice. To illustrate clearly the role 
of the phase diffusion in the interference pattern of the 
expanding subcondensates, we consider here an experi- 
ment scheme which can be realized in future experiment. 
Firstly, a cigar-shaped condensate is formed in a mag- 
netic trap with the harmonic angular frequencies ut x — 
2tt x 9 Hz and uj± — 2tt x 92 Hz at a temperature far be- 
low the critical temperature. The atoms are transferred 
into the lattice potential by increasing the power of the 
laser beam so that the harmonic angular frequency of the 
well is approximated to be 2tt x 14 kHz in a;— direction. 
For uj x — 2ir x 14 kHz, based on the well-known Bose- 
Hubbard model, the tunnelling time between neighboring 
subcondensates is of the order of 1 ms. In this situation, 
the array of subcondensates should be regarded as fully 
coherent and there is no independent random phase for 
the subcondensates in different wells. Although there is 
still a phase diffusion for the whole condensate, one can 
not observe this through the interference pattern. To ob- 
serve the effect of the phase diffusion, one can rapidly 
increase the lattice potential depth to a value of uo x = 
2tt x 52 kHz within a time (for example 50 /is) much 
smaller than the tunnelling time. In this situation, the 
phase of every subcondensate is the same after the lat- 
tice potential ramps to its final strength. For the lattice 
potential with ui x = 2tt x 52 kHz, the tunnelling time 
between neighboring wells is of the order of 1.5 s. Thus, 
for the holding time th of the optical lattice being much 
smaller than 1.5 s, the phase diffusion of the subconden- 
sates in different wells can be regarded as independent 
from each other. In this situation, there is a random 
phase for each subcondensate with the development of 



the holding time and this will lead to an obvious effect 
on the interference pattern after the combined potential 
is switched off. 

For the experimental parameters considered here, we 
have S 2 4> = 1-60 x th- Based on the result given by Eq. 
(20), the probability distribution of the phase for the 
subcondensate in each well can be obtained through Eq. 
(17). For the interference pattern to wash out completely, 
the time scale of the holding time is obtained by using 
5 2 <f> = ir 2 . In the experimental parameters used here, for 
the subcondensate in each well 5 2 <fi is still much smaller 
than 7r 2 even when the holding time th of the optical 
lattice is Is. Due to the fact that there are many subcon- 
densates induced by the optical lattice, however, we will 
show that there is a significant effect of the phase diffu- 
sion to the interference pattern after the optical lattice 
is switched off. 

To observe a notable effect of the phase diffusion, we 
consider the case when only the optical lattice is switched 
off after the holding time th = 0.2 s. Assume to is the 
time after the optical lattice is switched off. For the sub- 
condensate confined in the fcth well of the optical lattice, 
applying the propagator method used in Ref. [28] it is 
easy to get the normalized wave function iff. {x,to) after 
only the optical lattice is switched off: 



(fk{x,t ) = A n , 



x cxp 



sinw^o (ctgw^o + h) 



cxp[ii? fc (5<j>,t h )] 



(kdcosuj x ta — xY 



x exp 



2cr 2 sin 2 uj x t Q (ctg 2 a> x i + 7 2 ) 



i (kd cos u x to — x) ctgu^to 
27CT 2 sin 2 u> x t (ctg 2 uj x t + 7 2 



x cxp 



i [x 2 cosLd x t + k 2 d 2 cosoj x t — 2xkd) 
27CT 2 smtd x to 



(21) 



where A n = l/ir 1 / 4 /! 1 / 2 is a normalization constant, 
the dimcnsionless parameter 7 = h/muj x a 2 with a the 
width of the subcondensate in each well [28] . The factor 
Rk (5(j),th) represents a random phase due to the phase 
diffusion of the subcondensate in the fcth well. In numeri- 
cal calculations, Rk (S<j), th) is generated according to the 
probability distribution function P(A<j>,th) (see (17)). 
Then we obtain the density distribution in x— direction 
as follows: 



n(x,t ) = 5 



k M 

E 

k— — k]\ 



k 2 
k 



A I 



— ) VkixM] 



(22) 



where S = 15N a uk%[/ (l6fc|^ — l). N a u is the total num- 
ber of atoms in the array of condensates and 2k m + 1 is 
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the total number of the subcondensates induced by the 
optical lattice. In the present work, by using the experi- 
mental result in [20], ku is chosen to be 100. 

It is obvious that the density distribution n(x,to) 
at x = reaches the maximum value at time ti — 
(21 — l)ir/2u x , where I is a positive integer. Displayed 
in Fig. 1 is n(x,to) (with N a i{EA\ as a unit) at ti with 
the holding time th = 0.2 s. We see that the phase 
diffusion of the array of the subcondensates makes the 
central density of the interference pattern decrease sig- 
nificantly. Displayed in Fig. 2 is the density distribution 
n (x — 0, to =ti) versus the holding time th- A exponen- 
tial damping of the central density is clearly shown in the 
figure. 




-600 -400 -200 200 400 600 

x (unit of d) 

FIG. 1. Displayed is the density distribution of the inter- 
ference pattern at to = U for the case of the holding time th 
of the optical lattices being 0.2 s. The density distribution 
n(x,ti) is in unit of NaiiSA^, while the coordinate x is in 
unit of d. In this figure, the position of the cross x indicates 
the density at x = in the case of th = s. The figure shows 
clearly that the phase diffusion has the effect of decreasing 
significantly the maximum density of the central peak. 




It is straightforward to calculate the interference pat- 
tern of the expanding array of BECs when the magnetic 
trap and optical lattice are both switched off. Shown in 
Fig. 3 is the interference pattern for the holding time 
th = 0.2 s after the combined potential is switched off for 
30 ms. In comparison with the case of th = s (dashed 
line in Fig. 4), the strong noise in Fig. 3 shows that the 
phase diffusion has an important effect on the interfer- 
ence pattern of the array of subcondensates. The solid 
line in Fig. 4 displays the average density distribution of 
twenty ensembles for the holding time th = 0.2 s. We see 
that the maximum density of the central and side peaks 
decreases significantly due to the presence of the phase 
diffusion. Generally speaking, there are four parameters 
ui±, uj x , ku, th which are related closely to the effect of 
the phase diffusion. The present work shows that increas- 
ing these parameters will contribute to the observation 
of the phase diffusion for BEC in optical lattices. 
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x (unit of d) 

FIG. 3. After the holding time 0.2 s of the optical lattices, 
when both the magnetic trap and optical lattices are switched 
off, the density distribution of the expanding condensates is 
shown for to = 30 ms. In the figure, n(x,ti) is in unit of 
NaiiS, A^, while x is in unit of d. Compared with the density 
distribution of the zero holding time (dashed line in Fig. 4), 
the effect of the phase diffusion on the interference pattern is 
very obvious. 



FIG. 2. Shown is n(x = 0,ti,t h ) versus the holding time 
t h . The exponential damping of the central peak is clearly 
shown in the figure. 
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of China under grant Nos. 10205011 and 10274021. 
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FIG. 4. The solid line displays the average density distri- 
bution of twenty ensembles when both the magnetic trap and 
optical lattices are switched off for to = 30 ms in the case 
of the holding time th = 0.2 s. The dashed line shows the 
interference pattern in the case of the holding time th = 0. 



IV. DISCUSSION AND SUMMARY 

Based on the analysis of the collective excitations, we 
have investigated the phase diffusion at low tempera- 
ture in an array of subcondensates formed in a combined 
potential consisting of harmonically magnetic trap and 
ID optical lattice. We have shown that the effect of 
phase diffusion on the interference pattern is evident and 
thus can be measured from the expanding subconden- 
sates after the optical lattice is switched off. We point 
out that the maximum density of the interference pattern 
decreases significantly due to the mixing of the phase dif- 
fusion appearing in different subcondensates formed by 
the optical lattice. This suggests a way to detect experi- 
mentally the notable phase diffusion of BEC. For the case 
of BEC in 3D optical lattice, however, it seems that the 
theory developed here can not be used due to the fact 
that there are only an average atom number of up to 2.5 
atoms per lattice site [19]. 

Note that a recent experiment has been conducted in 
Rcf. [21] where the interference pattern is measured when 
only the optical lattice is switched off. Clearly this type 
of experiment can be used to test the theoretical predic- 
tions provided in this paper. We stress that for a single 
subcondensate the phase fluctuations increase very slowly 
with the development of the time. However, the phase 
diffusion plays an important role in the interference pat- 
tern of the expanding subcondensates, as shown in the 
present work. 
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